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Correlation effects on the static structure factor of a Bose gas 
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\l ' A theoretical treatment of the static structure factor S{k) of a Bose gas is attempted. The low 

^^ , order expansion theory is implemented for the construction of the two body density distribution, 

while various trial functions for the radial distribution function g{r) are used. g{r) introduces the 

atomic correlations and describes the departure from the noninteracting gas. The Bose gas is studied 

S!h ■ as inhomogeneous one, trapped in harmonic oscillator well, as well as homogeneous. A suitable 

3 ' parametrization of the various trial functions g{r) exists which leads to satisfactory reproduction of 

' , the experimental values of S(fc), both in inhomogeneous case as well as in homogeneous one. The 

O^ ' phonon range behavior of the calculated S{k) is also addressed and discussed both in finite and 

Cn ' infinite Bose gas. 
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C^ ■ I. INTRODUCTION 

S: 

1-^ ' Spectroscopic studies have been used to assemble a complete understanding of the structure of atoms and simple 

r^ . molecules [1]. The static structure factor S{k) is a fundamental quantity, which is connected with the atomic structure, 

Q ' and it is the Fourier transform of the radial distribution function g{r). S{k) gives the magnitude of the density 

O fluctuation in the system (atomic, molecular, electronic or nuclear system) at wavelength 27r/fc, where k is the 

momentum transfer. In recent papers, the Bragg spectroscopic method was used to measure S{k) either in the 

,__! [ phonon regime [1] or/and in the single-particle regime [2] . 

J> . More specifically, the character of excitations in a weakly interacting Bose-Einstein condensed gas depends on the 

^^ ' relation between the wave vector of the excitation k and the inverse healing length ^^^ — \/2mCsl'h = y^Sttop, which 






is the wave vector related to the speed of Bogoliubov sound Cs — yji/m, where /i — Anh ap/m is the chemical 
potential, a is the scattering length, p is the condensate density, and m is the atomic mass [1,3]. There are two 
^-^ . different regimes of excitation described by the Bogoliubov theory for the zero-temperature weakly interacting Bose- 
tH- ' Einstein condensate: the phonon and the free-particle regimes. For small wave vectors {k ^ £,~^), the gas responds 
f^ , collectively and density perturbations propagate as phonons at the speed of Bogoliubov sound. For large wave vectors 
^^ ■ {k ^ ^~^) the excitations are particle- like with a quadratic dispersion relation [1,3]. 

C^ [ In the present letter we focus on the theoretical calculation of S{k) of a Bose gas. All the calculations are performed 

H ■ in the zero-temperature limit (T = 0). Various types of trial functions g{r) are used and the entailed S{k) are compared 
with the experimental data. The correlation parameters of g(r) are determined by fit to the experimental data. The 
theoretical investigation is performed by applying two approaches. In the first approach, the Bose gas is treated as 
inhomogeneous gas, which is exactly the case of the experiments [2,1]. We consider that the Bose gas is trapped in 
Q , a harmonic oscillator well with condensate radius R. In the second approach, the trapped bose gas is treated, in a 
rough approximation, as homogeneous gas in regard to its elementary excitations. This approximate treatment is 
reasonable in the sense that the coherence length ^ satisfies the condition ^ <C i? [4] . Actually in the experiments of 
Refs. [1,2] this condition is well satisfied. 
5— i ' As mentioned before, the key quantity of the above approaches, is the radial distribution function g(r), which 

describes the relative probability for finding two particles at a distance r apart and consequently it introduces the 
atomic correlations. In general, the function g{r) can be calculated by the variational method [5]. However, in dilute 
Bose gas, in the framework of the lowest-order cluster expansion, trial expressions for g{r) can be used [5-7] and the 
parameters can be fitted to reproduce the experimental data of S{k). Up to now, S{k) in a trapped Bose gas, is 
predicted successfully (compared to the experimental data) by the local density approximation [8,2]. 

The motivation of the present work is the theoretical study of S{k), considering the Bose gas as a many-body 
system. It is well known that the ideal Bose gas model (noninteracting model) fails to reproduce the experimental 
values of S{k). We focus our treatment on the effect of interatomic interaction on the properties of S{k), both in 
small and large values of the momentum transfer k. It is fundamental to investigate how these effects modify the 
picture of S{k) of an ideal gas and to predict values consistent with the experimental data. The possibility of linear 
behavior of S{k) for small k, both in inhomogeneous and homogeneous Bose gas has also been examined. In section 
2 we treat the inhomogeneous Bose gas while the case of homogeneous Bose gas is investigated in section 3. 
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II. INHOMOGENEOUS BOSE GAS 



A dilute inhomogeneous Bose gas can be studied using the low-order approximation (LOA) [5-7]. In the LOA the 
two-body density distribution (TBDD) has the form 



p(ri,r2) = Cp{ri)p{r2).f{ri2) = Cp(ri)p(r2)g(ri2), 



(1) 



where /(ri2) {ri2 = r = |ri — r2|) is the Jastrow correlation function [6], g{ri2) is the radial distribution function, C is 
the normalization factor which ensures that / p(ri, r2)(iri(ir2 = N{N — 1) (N is the number of the atoms of the Bose 
condensate) and p(r) is the density distribution (DD) of the system. The TBDD is proportional to the probability 
of simultaneously finding an atom at ri and another at r2. In the uncorrelated case (noninteracting gas) the TBDD 
takes the simple form 



A^- 1 
p(ri,r2) = -^^p{ri)p{r2). 



(2) 



In the present work we consider that the atoms are confined in a harmonic oscillator trap where the ground state 
single-particle wave function V'o(r) has the form 



i'oir) 



7^3/453/2 



exp[-rV(26^)], b=[h/{muj)]^^'^. 



The normalization of the DD is / | p{r) p rfr = / | ipo(r) p dr = N. 
The static structure factor S{k), in a finite system, is defined as [i 



or using Eq. (1) 



S{k) = 1 + ^ / e*k(ri-r.) [p(i.i, r2) - p(ri)p(r2)] dndr2 
S{k) = 1 + 1 y e*'^(--i-'-^)p(ri)p(r2)[C5(ri2) - l]dr,dr2 



(3) 



(4) 



The integration in Eq. (4) can be performed if the function g(r) is know. g{r) must obey the rules g{r = 0) = and 
limr^oo g{f) —^ 1- The first rule introduces the repulsive correlations between the atoms and the second the absence 
of such correlations in long distances. In general the form of S{k) is affected appreciably from the form of g{r). More 
specifically, the long range behavior of g{r) affects S{k) for small values of k while the short range behaviour of it 
affects S{k) for large values of fc as a direct consequence of the Fourier transform theory. 

In the present work, we choose two trial forms for g{r). The first one is a gaussian type which has been extensively 
and successfully used for the study of similar problems in atomic physics (Bose gas, liquid helium) as well in nuclear 
physics. The relevant g{r) and the entailed S{k) are 



Case 1 



g{r) = 1 - exp[-/3r2] 

5(fc) = 1 + NiCi - 1) exp[-fcg/2] - 



NCi 



(l + 2y2)3/2 



exp[-fc2/2(l + 2y)]. 



(5) 



where kb = kb, y = f3b'^, (3 is the correlation parameter and Ci is the normalization factor. 
The second trial function g{r) and the relevant S{k) are of the form 



Case 2 
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S{k) = 1 + N{C2 - 1) exp[-A:bV2] 

NC2 
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2i5/2a67r/ch ^ 



i=l 



Aexph/3f/4] + V^(l 



Y)erf(-) 



(6) 



where a is the correlation parameter, a^ are known coefficients, f3i = f3i(kb,a,b), erf(z) = 2 
the normalization factor. 



dt and C2 is 



Trial functions g{r) of hard-sphere form {g{r) — l^a/r) and soft repulsive form {g{r) = 1 — sinar/ar) are also used. 
Both of them, lead either to the ideal Bose gas results (negligible correlations) or to abnormal fluctuating negative 
and positive values for S{k) (strong correlations) and consequently in complete disagreement with the experimental 
data [9]. Therefore, it should be emphasized that the choice of a trial function g{r) should be made under certain 
restrictions in order to get reasonable values for S{k). 

The behavior of S(k), in case 1, for different values of the correlation parameter f3 is shown in Fig. 1(a). It is 
obvious that the effect of correlations, induced by the function g(r), becomes large when the parameter /3 becomes 
small and vice versa. The case where /3 —^ oo, corresponds to the uncorrelated case (HO). For the values of k employed 
in the experiment of Ref. [2] (hereafter EXP-1) the prediction of the HO model is always close to 1 for S(k). When 
the correlation parameter (3 decreases considerably (strong correlations) the theoretical prediction of S{k) is in good 
agreement with the experimental data. The value /3 — 5.3 ^m'^^ gives the best least squares fit in that case. In general 
the gaussian form of g(r), in spite of its simplicity, reproduces fairly well the experimental data of EXP-1, both in low 
and high values of the momentum k. It reproduces also the experimental data of Ref. [1] (hereafter EXP-2) as can 
been seen from Fig. 1(c). Within our theoretical model, the gaussian type of g{r) is flexible enough to obtain values 
for S{k) in agreement with the experimental data. 

Fig. 1(b) displays the results in case 2, which are compared with those of the data of EXP-1. The model reproduces 
well the experimental data in the range 1.5 — 3 fj,m~^ (with best least squares fit value a — 1.34 ^m~^), but fails in 
the range fc > 3 nm~^. The main drawback of this model is the predicted negative values of S{k) in the range close 
to fc = when the correlation parameter a decreases considerably (strong correlation case) . 

The correlation function g{r) corresponding to cases 1 and 2 for the correlation parameters /3 — 5.3 fim~^ and 
a ~ 1.34 jim^^ respectively is sketched in Fig. 2(a). Those values of the parameters (3 and a give the best x^ value in 
the fit of the theoretical expressions of S{k) to the data of EXP-1. The most striking feature in case 2 is the existence 
of strong correlations, introduced by (?(r), in order to reproduce the experimental data of S{k). It is worthwhile to 
point out that g{r), in case 2, exhibits fluctuations in the range r > 2 /i?n but this is not visible in Fig. 2(a). 

The possibility of a linear dependence of S{k) on k for small values of k, as predicted from other works [8], is 
prohibitive, on the basis of Eq. (3) at least in the case where the trap is an harmonic oscillator one. That can be 
seen considering the ground state wave function to be the harmonic oscillator one and transforming ri and r2 in Eq. 
(3) into the coordinates of the relative motion (r = ri — r2) and the center of mass motion (R = (ri -f- r2)/2). After 
some algebra S'(k) takes the form 

S{\^)^ j ^'^^e-^\C9{r)-l]dv (7) 

For finite systems, as is a trapped Bose gas, we can expand the exponential e*'"', since r is bounded. So that: 

e*'^"- = l + 2kr+^^ + ^^ + --- (8) 

Substituting Eq. (8) into Eq. (7) and considering that the terms with odd powers of k do not contribute on the 
integral, S{k) takes the form 

5(fc)-aifc2+a2fc'' + --- (9) 

Thus, for small values of fc, S{k) depends linearly on fc^. The gaussian factor e"*" , originating from the harmonic 
oscillator wave function of the trapped Bose gas, ensures the convergence of the integrals a^ corresponding to the even 
powers of the expansion. 

III. HOMOGENEOUS BOSE GAS 

The condensate of an inhomogeneous (finite) Bose gas, can be treated as homogeneous in regard to its elementary 
excitations considering that the coherence length ^ satisfies ^ -^ R [4]. Actually, this is a rough approximation, but 
can effectively describe fairly well the excitation properties of a trapped Bose gas. 

In infinite systems p(ri) is constant (p(ri) = p) and thus the TBDD is given by p(ri,r2) ~ p^g{ri2). Thus, the 
structure factor of homogeneous gas (or quantum fiuid in general) is given by the relation [10] 

S(k) = l + p fe^'^'^[g{r)-l]dr, (10) 



or after performing the angular integration 

S{k) = 1 + 47rp / ^^^r^igir) - l]dr. (11) 

Jo ^^ 

The necessary conditions which must be obeyed by g{r) and S{k), when a given g(r) is not generated directly and 
explicitly by a wave function, include [10] 

Tik 
g{r) > 0, S{k) > 0, lim S{k) 



k^o 2mc 

f-OO 

S{Q) = ^ Anp I r^[g{r) - l]dr = -1 (12) 

JO 

The third condition of (12) imposes the condition for the phonon like excitation on the Bose gas and is satisfied if 
and only if [10] 

limr4[;^-l] = -^^ (13) 

where the bar denotes an average over a range Sr somewhat larger than p^^'^ [10]. 

We have used three different trial expressions for g{r) to study S{k) in a uniform Bose gas. The first two were 
used already in the study of the inhomogeneous Bose gas. However, the relevant structure factor, as expected, has 
different forms compared to inhomogenous case. The used forms of g{r) and the entailed S{k) are 

Case 3 

g{r) = 1 - exp[-/3r^] 

5(fc) = l-exp[-^] (14) 



Case 4 
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k <2a 

S{k) = '{ 2-f -A, 2a<k<4a (15) 

Aa < k 



Ha ' 
2a 



Case 5 is a combination, in a way, of cases 3 and 4 and takes the form 

Case 5 

sin ar 

g{r) ~ L exp[— crj 

ar 

S(k) = fc^(fc^ + 2c^-2a^) 

^' (a2 + c2)2 + fc2(p + 2c2-2a2) "^ "' 

It is worth noticing that in homogeneous cases, due to the fourth condition of Eq. (12), the correlation parameters 
oi g{r), in all the cases, depend directly on the density p. Thus, the fourth condition of Eq. (12), leads to the relations 
p = /33/2/7r3/2 (case 3), p = a^/n^ (case 4) and p ^ {a'^ + c'^)^/8cn (case 5). 

In the present work we compare also our results for S{k) with those of the hard-sphere interaction in uniform dilute 
gas, which was predicted long ago by Lee et al. [11] and has the form 

5(fc) = \ (17) 

y' k'' + Ibirap 

where a is the hard-sphere diameter (see also Refs. [12,13]). 

The values of S{k) in case 3 (gaussian g{r)) for various values of the constant density p (or of the correlation 
parameter /3) are displayed in Fig. 3(a). As in the corresponding case of the inhomogeneous Bose gas that type oi g{r) 
reproduces the data of both experiments, see also Fig. 1(c), when the density of the gas is approximately « 1 pm~^ — 



10^^ cm^^. Wc note that in the uncorrelated case (nonintcracting model), the static structure factor of homogeneous 
Bose gas has the constant value S{k) = 1. When the density of the Bose gas increases, the calculated values of 
S{k) deviates from the experimental data, exhibiting considerably lower values at the same range of the momentum 
transfer k. Moreover, S{k) behaves quadratically in the phonon regime (linifc^o S{k) ~ fc^), in contradiction with the 
condition (12), as a result of the improper long range behavior oi g{r) (see Eq. (13)). 

The behavior of S{k), in case 4, is illustrated in Fig. 3(b). A striking future in this case is the linear behavior of S{k) 
(limfe_,o S{k) ^ fc), in the phonon regime, as a consequence of the proper long range behavior of the function g(r). The 
best fit value of the correlation parameter a corresponds, as in the previous case, to density p ^ 1 p.m~'^ ~ 10^^ cm~'^. 
Increasing the density, lower values of S{k) appear for the same range of the momentum k. 

S{k) in case 5, is plotted in Fig. 3(c). In this case g{r), combines in a way cases 3 and 4. S{k) behaves quadratically 
in phonon regime {S{k ^ 0) ~ fc^) but reproduces quite well the experimental data, especially in the range 0—4 fj,m~^. 
S{k), derived by Lee et al. [11], is also plotted in Fig. 3(c). In that case, the phonon behavior is ensured, while the 
prediction is in good agreement with the experimental data. 

The theoretical results of S{k) in cases 1,3,4 and 5 using the best fit values of the parameters are compared with the 
experimental data of Exp-2 in Fig. 1(c). It is seen that in all cases there is a good agreement between the theoretical 
values and the experimental data. 

The behavior of the correlation function g{r), for cases 3,4 and 5, is sketched in Fig. 2(b). The most striking feature 
is that though the behavior of g{r) is almost the same in cases 3 and 4, the corresponding S{k), display different 
behavior, especially for the lower values of the momentum k. It is concluded that S{k) is very sensitive to the form of 
g{r) and specifically, the long range behavior of g{r) affects considerably the behavior of S{k) in the phonon regime 
(small values of k). This is a well known property of the Fourier transform. 

It is worth remarking that a correlation function of the hard-sphere form, g(r) = 1 — a/r, thought it is as realistic 
one, mainly for low values of the interatomic distance r, leads to S{k) = 1 — Airpa/k"^ which diverges when fc ^ 0. 

In conclusion, we report a theoretical calculation of the static structure factor S{k) both for inhomogeneous and 
homogeneous Bose gas, in the framework of the low order expansion theory, by applying various trial forms for the 
radial distribution function g{r). We compared our results with recent experimental data concerning trapped Bose 
gas. The correlation parameters of g{r) are adjusted in order to reproduce the experimental data. By applying 
suitable parametrization the experimental data are reproduced quite well. The low k behavior of the calculated S{k) 
is also addressed and discussed both in inhomogeneous and homogeneous Bose gas. 
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correlation function, N. Katz for literature information about the experimental data for S{k) and Prof. A. Herrera- 
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FIG. 1. The static structure factor S'(fc) of the inhomogeneous and homogeneous Bose gas in various cases versus the 
momentum fc, (a) in Case 1 for various values of the correlation parameter /3 as well as for the uncorrelated case (harmonic 
oscillator), (b) in Case 2 for the least squares best fit value of the parameter a, (c) in Cases 1,3,4 and 5 for the best fit values 
of the correlation parameters. The experimental points EXP-1 and EXP-2 are from references [2] and [1] respectively. For the 
various cases see text. 
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FIG. 2. The radial distribution function gir) for Case 1 and 2 (a) (corresponding to inhomogeneous Bose gas), and Cases 
3,4 and 5 (b) (corresponding to homogeneous Bose gas) with the best fit values of the correlation parameters. 
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FIG. 3. The static structure factor S(fc) for Case 3 (a) and Case 4 (b), for various values of the density p. The Case 5 (c) for 
the best fit value of the density and also for the hard-sphere interaction of Lee et.al. [11] (indicated LHY). The experimental 
data are from Ref. [2] 



